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We present the first-principles studies on the hydrostatic pressure effect of the electronic properties 
of black phosphorus. We show that the energy bands crossover around the critical pressure P c — 1.23 
GPa; with increasing pressure, the band reversal occurs at the Z point and evolves into 4 twofold- 
degenerate Dirac cones around the Z point, suggesting that pressured black phosphorus is a 3D Dirac 
semimetal. With further increasing pressure the Dirac cones in the T-Z line move toward the T point 
and evolve into 2 hole-type Fermi pockets, and those in the Z-M lines move toward the M point 
and evolve into 2 tiny electron-type Fermi pockets, and a band above the Z-M line sinks below E f 
and contributes 4 electron-type pockets. A clear Lifshitz transition occurs at P c from semiconductor 
to 3D Dirac semimetal. Such a 3D Dirac semimetal is protected by the nonsymmorphic space 
symmetry of bulk black phosphorus. These suggest the bright perspective of black phosphorus for 
optoelectronic and electronic devices due to its easy modulation by pressure. 


I. MOTIVATIONS 


The van der Waals force between layers makes layered 
compounds graphite, M 0 S 2 , WTe 2 and TaS 2 , etc.,be eas¬ 
ily exfoliated to monolayers as two-dimensional materi¬ 
als, such as graphite to graphene and M 0 S 2 to its mono- 
layer. As a possible candidate of optoelectronic and elec¬ 
tronic material, the energy gap of black phosphorus varies 
from 0.3 eV in bulk to about 2 eV in monolayer, filling 
the optical interval between small energy gap 0-0.3 eV in 
graphene and large energy gap of 1-2 eV in semiconduc- 
tive dichalcogenides [1]. Recently, black phosphorus in 
bulk, multilayers and monolayers have received consider¬ 
able interests [2-4]. Zhang et al. demonstrated that mul¬ 
tilayer black phosphorus exhibits as high as 10 5 drain cur¬ 
rent modulation and 10 3 cm 2 V~ 1 s~ 1 charge mobility (U, 
showing that black phosphorus thin film might be a good 
potential candidate for field effect transistor [6]. More re¬ 
cently it is shown that a moderate hydrostatic pressure 
could not only drive semiconductive black phosphorus 
to metallic, but also tune multiple Fermi surfaces and 
Lifshitz point in the magnetotransport and Shubnikov- 
De Haas oscillation measurements [7|. These properties 
make black phosphorus a new exciting field both in ma¬ 
terial sciences and condensed matter physics j8|. 

Hydrostatic pressure could easily change the crys¬ 
tal structure of black phosphorus and modify its elec¬ 
tronic properties 1, 0. As the most stable allotrope 
of phosphorus, black phosphorus exhibits three differ¬ 
ent phases under moderate high hydrostatic pressures: 
the orthorhombic phase with wrinkled hexagons, the 
graphene-like rhombohedral phase with hexagonal lattice 
for the pressure P>4.5 GPa 0, and the simple cubic 
phase for P>10.3 GPa 0. Under high pressure the 


P-P bond lengths and bond angle of black phosphorus 
display strong pressure dependence, and exhibit highly 
anisotropy in compressibility [9]. Though a strain may 
crucially modify the energy gap of monolayer black phos¬ 
phorus |3[ , it is completely different from the hydrostatic 
pressure effect on bulk black phosphorus, since the van 
der Waals force leads to considerable interlayer coupling 
a One anticipates that a high hydrostatic pressure sig¬ 
nificantly compresses the interlayer distance of P atoms, 
however, it remains a puzzle and unclear scenario: how 
do the electronic structures evolve with increasing hydro¬ 
static pressures? Especially, how does the Lifshitz point 
occur under pressure? 


On the other hand, searching for 3D Dirac semimet¬ 
als has attracted great attention in recent years 12 M■ 
Recently Kim et al claimed 0 that they observed 
the Dirac semimetal state by depositing K ions on the 
surface of black phosphorus, however, their Dirac cones 
near the T point actually always sink below Ep. Cava 
et al. [13] pointed out that a crystal with glide planes 
and scream axes might form a 3D Dirac semimetal when 
the electronic structures of the crystal possess a twofold- 
degenerate bands, since the nonsymmorphic symmetry 
in the crystal lattice leads to the band inversion at a 
point in the Brillouin zone 00. The space group of 
bulk black phosphorus has a glide plane, thus it is a po¬ 
tential candidate of a 3D Dirac semimetal. To elucidate 
electronic and optical properties of black phosphorus un¬ 
der hydrostatic pressure, in this Letter, using the first- 
principles electronic structures calculations and analytic 
studies, we for the first time present the evolutions of 
the band structures and Fermi surfaces of black phospho¬ 
rus under hydrostatic pressure, and clearly demonstrate 
the Lifshitz transition from a semiconductor to a Dirac 
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semimetal. We argue that the non-symmorphic symme¬ 
try in crystal structure of black phosphorus protects it 
as a 3D Dirac semimetal. 


II. COMPUTATIONAL STRUCTURES AND 
PHONON SPECTRA 

II. 1 Computational details 


To explore the evolutions of crystal structures, elec¬ 
tronic structures and physical properties with hydrostatic 
pressures, we first achieve the stable crystal structures 
of black phosphorus under various pressures within the 
framework of density functional theory (DFT) using the 
projector augmented wave (PAW) method [17!, Il3j as ini- 
plemented in the Vienna Ab initio Simulation Package 
(VASP) code 0. These structures are optimized by us¬ 
ing the optB88-vdW method 2^, El], so as to take into 
account the interlayer van der Waals force in bulk black 
phosphorus. 

The experimental lattice parameters of bulk black 
phosphorus are taken as the initial structure. In our 
work, the directions of lattice constants a and b are set to 
be along the zigzag (the y axis) and the armchair ( the 
x axis) directions, respectively, in the crystal structure 
of bulk black phosphorus; and c is along the interlayer 
direction (the z axis), as displayed in Fig. [lj Both the 
shape and volume of each supercell have been relaxed 
fully at a target pressure and all atoms in the supercell 
were allowed to move until the residual force per atom 
was drop below 0.001 eV/A. The energy cutoff for the 
plane-wave basis was set to 500 eV for all calculations. 
A k-mesh of 10 x 4 x 8 in band structure calculation was 
adopted to sample the first Brillouin zone of the primi¬ 
tive unit cell of bulk black phosphorus. Uniform k-mesh 
of 30x30 x 30 was adopted to obtain Fermi surfaces when 
black phosphorus changes into semimetal. Phonon cal¬ 
culations were carried out using the finite displacement 
method and PHONONPY package 
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|. A 3x3x4 

supercell and a 5x5x5 k-mesh were applied to calculate 
the force constants. 



FIG. 1: Crystal structure of bulk black phosphorus (a) and 
the first Brillouin zone and some high symmetric points of 
bulk black phosphorus (b). 


II.2 Crystal structure evolution and phonon spectra 


We find that the calculated lattice constants of bulk 
black phosphorus in the orthorhombic phase under var¬ 
ious pressures shown in Fig. [2] are in agreement with 
the experimental results 24|, |25j very well. In detail, the 


calculated lattice constants a and b are very close to the 
experimental results [24. 25], though the lattice constant 
c is slightly larger than the measured within 0.2%, since 
the GGA functionals with the udW-correction overesti¬ 
mates the layer distance c, implying that our theoreti¬ 
cal structural parameters are highly reliable. Further¬ 
more, the phonon spectra for our theoretical structures 
also demonstrate that the pressured structures are sta¬ 
ble. Our numerical results have shown that all of our 
theoretical structures have not any imaginary frequency, 
as shown in Fig. lAll in the APPENDIX. 


III. ELECTRONIC STRUCTURES EVOLUTION 
UNDER PRESSURE 


Electronic band structures were obtained by the VASP 
method based on the different pressurized structures. All 
of the calculations are cross-checked by using the Quan¬ 
tum Espresso [26| and WIEN2K codes, respectively, 
and the results are consistent with each other. Within 
the hybrid functional of Heyd, Scuseria, and Ernzerhof 
(HSE) framework we obtain an energy gap of 

0.34 eV at ambient condition, in agreement with the ex¬ 
perimental and recent theoretical data EE!. 

To save the computation resource and accord with the 
experiments 0, throughout this paper we present the nu¬ 
merical results within the modified Becke-Johnson (mBJ) 
exchange-correlation potential 0E1 for the semicon- 
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FIG. 2: Lattice constants of bulk black phosphorus under 
various pressures. The experimental data from Iwasaki et al 
[24] and Cartz et al. [25] are also shown as open triangles and 
squares for comparison. 
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ducting phase, and the PBE exchange-correlation poten¬ 
tial [32| for the metallic phase. It is now recognized that 
similar to the HSE or GW methods, the mBJ method 
could well describe a wide range of materials properties, 
including the band gap of semiconductors, while it some¬ 
times make false predictions of fundamental properties of 
metals [33]. During the application of hydrostatic pres¬ 
sure, the black phosphorus undergoes from semiconduc- 
tive phase to semimetallic one. It is desirable and self- 
consistent to adopt the same first-principles method with 
the same exchange-correlation potential to calculate the 
electronic properties of black phosphorus over wide pres¬ 
sure range. However, due to the limit of each method, 
we have to adopt different exchange-correlation poten¬ 
tial to calculate the energy band of the semiconductor 
and semimetal. 

The pressure dependence of band structures for the 
semiconductive (mBJ results) and for the semimetallic 
(PBE result) phases of bulk black phosphorus is shown 
in Fig. [3l It shows that the system opens a gap about 
0.28 eV at 0 GPa. With increasing the pressure to 0.5 
and 1.0 GPa, up to about P c =1.23 GPa, the band gap 
gradually closes and the two bands below and above 
E p touch at the Z-point. Thus the critical pressure for 
the semiconductor-semimetal transition is about P c =1.23 
GPa in the mBJ functional, which is comparable with 
the value observed in the recent experiment [7]. Since 
the energy bands begin to anti-inverse at P c , the black 
phosphorus becomes a semimetal with a degenerate Dirac 
point at the Z point, as seen the band structures in Fig. 

ETc). 

When the hydrostatic pressure increases up to 1.5 GPa 
(not shown), the energy bands already cross with each 
other around the Fermi energy Ei?, the anti-inversed 
energy bands develop the linear dispersions and Dirac 
points in the Z-M and Z-Y directions, respectively. We 
attribute to that p z orbitals along the c axis gradually 
overlap with increasing the pressure, as shown in Fig. 0J 
majorly contributing to the band anti-inversion. At a 
pressure of 2.0 GPa shown in Fig. (3][d), the Dirac cones 
in the Z-Y direction shift toward the Y point and de¬ 
velops two hole Fermi pockets. On the other hand, the 
Dirac cones in the Z-M line shift toward the M point and 
above E^ in energy, contributing 2 tiny electron pockets; 
meanwhile a band above the Z-M line sinks below Ef 
and develops 4 electron Fermi pockets. At pressures up 
to 2.5 ^ 4 GPa we studied, the volumes of the hole Fermi 
pocket gradually expands, and an electron Fermi pocket 
in the Z-M line becomes considerably large, in agree¬ 
ment with recent experimental observation p]|. Thus it 
is naturally expected that the hole and electron carriers 
should coexist under high pressure until 4 GPa, since high 
hydrostatic pressure does not break the charge balance. 

As well known, the wide bandwidth 2 p z orbitals of 
carbon atoms contribute the six Dirac cones of graphene. 
One may wonder which orbital or which a few orbitals 



FIG. 3: (Color Online) Evolution of energy band structures of 
bulk black phosphorus with increasing pressures from 0 GPa 
(a), 0.5 GPa (b), 1.23 GPa, the critical pressure for band 
crossover (c), 2.0 GPa (d). The first three figures and the last 
one are the mBJ and PBE results, respectively. 


contribute the Dirac cones when the hydrostatic pressure 
P>P C . After projecting the weight of orbitals, we could 
find that the 3 p z orbitals contribute the dominant weight 
of the Dirac cones. As seen in Fig. |4j under the high 
hydrostatic pressure, the 3 p z orbitals overlap with each 
other, and lead to the linearized spectrum. 



FIG. 4: (Color online) Bands structure with orbital projec¬ 
tion of bulk black phosphorus at 2.0 GPa (a) and schematic 
diagram of 3 p z orbital distribution contributed to the Dirac 
cones (b). 
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IV. DIRAC CONES AND 3D DIRAC 
SEMIMETAL: 

From the band structures one finds that when the 
hydrostatic pressure increases to the critical value of 
P c — 1.23 GPa, the valence band and the conduction band 
touch together at the high-symmetric Z point of the 
Brillouin zone and gradually intersect with each other. 
The energy band structures of the system exhibit anti¬ 
inversion under further pressure. As described above, 
when the pressure becomes larger than the critical point, 
P > P c , the anti-inverted points at the high-symmetric 
Z point evolve into Dirac points around the Z point, and 
this fourfold-degenerate point gradually splits into four 
Dirac cones with the increase of pressure, as sketched in 
Fig. [5] and confirmed by Fermi surfaces evolution under 
pressures shown later. Taking into account the spin de¬ 
gree of freedom, each Dirac cone remains twofold degen¬ 
erate. Note that these Dirac cones are twofold-degenerate 
once taking into account the spin freedom of degree, thus 
the influence of magnetic field deserves for further inves¬ 
tigation. 

The 3D Dirac semimetal can further be verified from 
the density of states (DOS) of black phosphorus near 
the critical pressure P c shown in Fig. [6l One ob¬ 
serves that the DOS vanishes near E p, suggesting well- 
defined semimetallic nature. These properties strongly 
indicate that bulk black phosphorus under intermediate 
pressure is the first realistic example of pressure-modified 
3D Dirac semimetal. At P— 3.0 GPa shown in Fig. [6j fi¬ 
nite DOS appears near Ep due to the upward and down¬ 
ward of bands, implying that the system becomes metal¬ 
lic. However, the DOS near Ep does not increase consid¬ 
erably. 

The appearance of the band inversion, Dirac points 
and Dirac cones in black phosphorus are particularly in¬ 
teresting, since in this calculation, neither considerable 
spin-orbit coupling is taken into account, nor a graphene- 
like honeycomb structure is realized in black phosphorus , 
though the crystal structure of black phosphorus changes 



FIG. 5: (Color Online) Sketched distribution of 4 Dirac cones 
in the first Brillouin zone of bulk black phosphorus under the 
pressure of 2.0 GPa. The two Dirac cones lie on the upper 
and lower surfaces of the Brillouin zone are identical. 



FIG. 6: (Color Online) Density of states (DOS) of bulk black 
phosphorus under the pressure of P c — 1.23 GPa (blue solid 
line) and 3.00 GPa (red dashed line). 

significantly. This raises a question: what is the physical 
origin of the formation of the Dirac semimetal? Consid¬ 
ering the space-group symmetry operators of the black 
phosphorus, we attribute the origin of Dirac points to 
the symmetry of the nonsymmorphical space group in 
black phosphorus. In the hydrostatic pressure range we 
studied, the space group of the crystal structure of black 
phosphorus is always Cmca, which possesses a nonsym¬ 
morphical operation of a glide along the c-axis. Under 
these symmetrical operations U, the Hamiltonian of the 
black phosphorus, H(k), satisfies: 

H( k) = U*(k)H(Uk)U(k). (1) 

At k z = 7r/c, or the Brillouin zone surface containing 
the Z point (0, 0, 7r/c), the energy band structures of a 
nonsymmorphical crystal may be connected, leading to 
the degeneracy Q- Such a connection and degeneracy 
lead to the fourfold-degenerate Dirac points [131, |35[ . 


V. FERMI SURFACES EVOLUTION WITH 
PRESSURES: 

From the preceding studies we have shown that crit¬ 
ical pressure P c is a Lifshitz point, which the topology 
of Fermi surface of bulk black phosphorus changes cru¬ 
cially. The Fermi surfaces and its cross-section of black 
phosphorus under the pressure of 2.0 GPa are displayed 
in Fig. [71 There are three different kinds of Fermi pock¬ 
ets: two hole pockets arising from the upward shift of 
two Dirac cones in the T-Z line above E_p, and two tiny 
electron pockets on the upper and lower surface of the 
Brillouin zone arising from the sink of two Dirac cones 
in the M-Z line below E p; another four electron-type 
pockets in the Z-M line come from the sink of the bands 
below Ep. These four Fermi pockets do not originate 
from the Dirac cones. 

The cross-section Fermi surface on the plane (0, k y , 
k z ) at 2.0 GPa is displayed in Fig. [71(b) . This shows that 
the Dirac points and the centers of the Fermi pockets 
are on the k x =0 plane. The extremal cross-section area 
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TABLE I: Extremal cross-section area Sf of electron and hole pockets at pressure 2.0 GPa. /3, a and a* stand for a Dirac 
dispersion hole pocket at Z-Y line, a normal electron pocket and a Dirac dispersion electron pocket at Z-M line, respectively. 


Pockets type 

/3 (hole) 

a (electron) 

a * (electron) 


1.06 a 

0.4 a 

0.23 a 


(0.42) b 

(0.22) b 

- 


- 

- 

(0.16) c 

Inequivalent numbers 

2 

4 

2 

Total 

2.12 

1.60 

0.46 


a Ours, the plane (0, k y , k z ) 

b Ref. [7], the plane ( k x , k y , 0), with the magnetic field H =20 T parallelled to 
the c-axis and under the pressure of 2.0 GPa. 
c Ref. [7], the plane ( k x , k y: k z ), where the direction of k x is 18° away from the 
k x direction, with H=20 T and under the pressure of 2.4 GPa. 


TABLE II: Effective masses of black phosphorus under various hydrostatic pressures and comparison with Ref. mm- 


P(GPa) 

m*/mo 

o 

* 55 

m*/mo 

0.00 

0.14(0.12[11], 0.08[37]) 

1.26(1.15[11], 1.03[37]) 0.16(0.15[11], 0.13[37]) 

e 0.50 

0.09 

1.24 

0.14 

1.00 

0.05 

1.21 

0.13 

0.00 

0.12(0.11 [11], 0.07[37]) 

0.90(0.71[11], 0.65[37]) 

0.34(0.30[11], 0.28[37]) 

h 0.50 

0.08 

0.82 

0.32 

1.00 

0.04 

0.78 

0.31 


Sf of the electron and hole pockets in the plane (0, k y , 
k z ) are estimated. Our results show Sf are 0.40, 0.23 
and 1.06 nm ~ 2 for a large electron pocket (<a), a tiny 
electron pocket (cC) and a hole pocket (/3), respectively, 
listed in Table I. As for Shubnikov-De Haas oscillation 
experiment jzj, it reports that 0.22 for a and 0.42 nm ~ 2 
for /3 in the plane ( k x , k y , 0), with the magnetic field 
H = 20 T parallelled to the c-axis and under the pres¬ 
sure of 2.0 GPa; while Sf for a * is 0.16 nm ~ 2 in ( k x , k y , 
k z ) with H = 20 T and under the pressure of 2.4 GPa. 
Our theoretical results have the same order magnitude 
with the experiment 7]. With consideration of the num¬ 
ber of two types of pockets (2 for /3, 4 for a and 2 for 
a * ) in the first Brillouin zone, the total 5ps of electron 
pockets and hole pockets are approximately equal to 2.00 
nm ~ 2 . The numbers of two kinds of carriers keep iden¬ 
tical, indicating a complete compensation. These theo¬ 
retical data show that the system is semimetal state, in 
good agreement with the results by the field dependence 
of transverse magnetoresistance experiment [7]. Also, the 
Dirac fermions and normal fermions coexist in the black 
phosphorus under the pressure larger than P c , providing 
a complete novel quantum system. 

The pressures can modify the electronic structures, 
leading to the interesting evolution of Fermi surfaces, 
and driving significant change in the topology of elec¬ 
tronic structures, i. e., the Lifshitz transition . We 
present Fermi surfaces topology of bulk black phospho¬ 
rus under different pressures in Fig. [8l It is shown that 
Fermi surface is null before the critical pressure for the 
semiconductor-semimetal transition. With increasing the 
pressure to 1.5 and 3.0 GPa, up to 4.0 GPa, the volumes 



FIG. 7: (Color Online) Fermi surface (a) and its projection (b) 
of bulk black phosphorus under the pressure of 2.0 GPa. The 
projection on the plane (0, k y , k z ) is shown. The two pockets 
in the Y-Z line are hole type, and the others are electron type. 


of hole Fermi pockets in Y-Z line gradually expands, and 
electron Fermi pockets in Z-M line grow largely, in con¬ 
sistent with the resent experimental results [7[ . The vol¬ 
umes of electron and hole pockets become large with in¬ 
creasing pressure, indicating that more and more hole 
and electron carriers are produced in the system. 
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P (GPa) P (GPa) 

FIG. 9: (Color Online) Dependence of the average Fermi 
velocity of bulk black phosphorus under the pressures from 1.5 
GPa up to 4.0 GPa by the PBE (a) and the mBJ potentials 
(b), respectively. 


FIG. 8: Calculated Fermi surfaces for bulk black phosphorus 
from 1.23 GPa (a), 1.5 GPa (b), 3.0 GPa (c) and 4.0 GPa (d). 


VI. ANISOTROPIC PHYSICAL PROPERTIES 


The anisotropic electronic structures in black phos¬ 
phorus result in highly orientation-dependent properties, 
such as the effective masses of carriers and Fermi veloci¬ 
ties, etc. It displays the hydrostatic pressure dependence 
of the Fermi velocities (Vf) of the hole and electron car¬ 
riers by the PBE potential shown in Fig. [9^a) and the 
mBJ one in Fig. (9^b), respectively. We focus on the 
Fermi velocity of Dirac cones near the Fermi level, which 
is important for the fundamental properties of the sys¬ 
tem. Note that we calculate the average Fermi velocity 
from two branches of Dirac cone. It has shown that the 
considerable pressure dependence of the Fermi velocity 
both for the PBE and the mBJ results, while, the pres¬ 
sure dependence of Vf by the mBJ potential is signifi¬ 
cantly larger than that by the PBE potential. Of course, 
more elaborated experiments are expected to confirm our 
theoretical details in the future. 

In addition, It is seen that the Fermi velocities mono- 
tonically increase with increasing the hydrostatic pres¬ 
sure. The velocities along the Z-Y direction are almost 
twice of those along the Z-M direction, showing strong 
transport anisotropy. Meanwhile, the effective masses of 
the hole and electron carriers also display considerable 
anisotropy. As listed in Table II, with increasing hydro¬ 
static pressure, the effective masses of electron carriers 
in the y- and z-axis almost do not alter, and the effective 
masses of hole carriers along the z-axis change slightly. 
However, the effective masses along the x-axis decrease to 
about one-third of the ambient pressure, demonstrating 
strong pressure modulation. 


VII. SUMMARY 


This work presents the nontrivial influence of hy¬ 
drostatic pressure on bulk black phosphorus, which 
is distinct different from the strain-induced effect [3]. 
The hydrostatic presure shortens the distance of two P 
atomic layers and the height of P atomic layers, which 
results in considerable changes in electronic structures 
and Fermi surface topology, hence a Lifshitz transition. 
It also leads to the appearance of 4 twofold-degenerate 
Dirac cones at a critical pressure P c , implying that 
black phosphorus under hydrostatic pressure becomes a 
3D Dirac semimetal. More recent work suggested that 
the these Dirac points may form unusual nodal ring 
in pressured black phosphorus [36| These show that 
pressure controlled layered black phosphorus may have 
potential application in optoelectronic and electronic 
device. 


The authors L.-J. Zou and D.-Y. Liu thank X.-G. Gong 
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Physics, Fudan University. This work was supported by 
the NSF of China under Grant No. 11474287, 11274310, 
11104274 and the key project 11534010. Numerical cal¬ 
culations were performed at the Center for Computa¬ 
tional Science of CASHIPS. 


VIII. APPENDIX 


To make sure that pressured structures are thermody¬ 
namically stable, we calculate the phonon spectra of bulk 
black phosphorus for our theoretical structures under the 
hydrostatic pressures from 0 GPa to 4 GPa, as shown in 
Fig. IA11 Our results show that the phonon spectra of 
these structures have no any imaginary frequency, im¬ 
plying that the optimized geometries really locate at the 
minimum point of the potential surface. 
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FIG. Al: Calculated phonon dispersion of bulk black phos¬ 
phorus under pressures of P— 0, 1, 2, 3, and 4 GPa. 
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